Centrality Dependence of Hadron Multiplicities 
in Nuclear Collisions in the Dual Parton Model 



A. Capella and D. Sousa 

Laboratoire de Physique Theorique^ 
Universite de Paris-Sud, Bailment 210, F-91^05 Orsay Cedex, France 



Abstract 

We show that, even in purely soft processes, the hadronic multiplicity in nucleus- 
nucleus interactions contains a term that scales with the number of binary collisions. 
In the absence of shadowing corrections, this term dominates at mid rapidities and 
high energies. Shadowing corrections are calculated as a function of impact pa- 
rameter and the centrality dependence of mid-rapidity multiplicities is determined. 
The multiplicity per participant increases with centrality with a rate that increases 
between SPS and RHIC energies, in agreement with experiment. 
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1 Introduction 



Hadron multiplicities per unit rapidity at SPS energy show an approximate scaling 
with the number of participants. This property is known as "wounded nucleon model" 
(WNM) [|TJ. However, a precise determination of the ratio dN ch / dr](r]* = 0)/n part , both 
at CERN-SPS @] and at RHIC ||, shows a steady increase from peripheral to central 
collisions. This increase is larger at RHIC energies, where the data show no saturation 
for the most central bins. 

It is a widespread belief that "soft" processes lead to a scaling in n part while "hard" 
ones lead to a scaling with the number n of binary collisions. While, this is the case for 
total cross-sections, it is not so for single particle inclusive cross-sections. If one neglects 
the effects of shadowing (i.e. nuclear effects in structure functions associated to triple 
Pomeron interactions), there is a theorem known as Abramovski, Gribov, Kancheli (AGK) 
cancellation Q, based on general principles, according to which da ch /dy scales with n at 
mid-rapidities and asymptotic energies. The best way to understand this cancellation is 
to illustrate it with a model that satisfies these general principles, namely the Glauber 
model. 

Let us consider for simplicity pA scattering. The cross-section for n inelastic colli- 
sions of the proton with n nucleons of the target is given by the probabilistic expression 

< A (b) = ( A ) U PP T A (b)Y(l - a pp T A (b)) A ~ n . (1) 



where T A (b) are nuclear profile functions normalized to unity. Using this equation it is 
easy to see that 

<= U%j:ai A {b)^A^ . (2) 

J n=l 

However, for the non-diffractive inclusive cross-section one has 

dnP A r A 

^— oc d%Y, n < A ( b ) <* Al ■ ( 3 ) 
dy J n =l 

Thus, one obtains the behaviour A 1 typical of hard processes. 

The derivation of (Bp assumes that the inclusive cross-section for n inelastic col- 
lisions is n times the corresponding one for a single collision (or, more generally, that 
the hadronic plateau produced in an inelastic collision does not depend on the number 
of inelastic collisions). Clearly, such a condition is only true at mid-rapidities and for 
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asymptotic energies. A basic idea in the WNM is that a nucleon (here the projectile), 
when wounded once, looses its ability to produce extra particles in further collisions. This 
may be valid at low energies, when the projectile undergoes successive collisions with the 
nucleons of the target. However, at high energies, the space-time development of the 
interaction implies that the n collisions are "parallel", i.e. they originate from different 
constituents of the projectile wave function and take place simultaneously. (Technically, 
it means that one is dealing with non-planar diagrams). This is the so-called Glauber- 
Gribov model ||. A more detailed discussion can be found in |J. 

As discussed above, at finite energies the constraints of energy-momentum conser- 
vation lead to a violation of the AGK cancellation. The dual parton model (DPM) |jj 
and the quark gluon string model (QGSM) ||, while obeying AGK cancellation at mid 
rapidities and asymptotic energies, contain an "educated guess" on its violation at finite 
energies. These models are based on the quark-gluon content of hadrons in the framework 
of the Glauber-Gribov model, and on the large N expansion of non-perturbative QCD. 
The charged multiplicities per unit rapidity are given by [0, [| 



dN' 



ch 

AA, 



(y, b) = n A (b) N« "* (y) + iV* ~ qq (y) + (2k - 2)N« 



q s -q s 



+ 



dy 

(n(b)-n A (b))(2kN^(y)) . (4) 
Here P and T stand for the projectile and target nuclei, 



n(b) = a pp J d 2 s A 2 T A (s) T A (b - s)/a AA (b) = a pp A 2 T AA (b) /a AA (b) (5) 
is the average number of binary collisions and 



n A (b) = n part (b)/2 = J d 2 s A T A (s) [1 - exp(-cr pp A T A (b - s)} /a AA (b) (6) 

the average number of participants of nucleus A. k is the average number of inelastic 
collisions in pp and /i(b) = kn(b)/n A (b) is the average total number of collisions suffered 
by each nucleon. The first term in (Q) is the plateau height in a pp collision, resulting from 
the superposition of 2k strings, multiplied by n A . Since in DPM there are two strings per 
inelastic collision, the second term, consisting of strings stretched between sea quarks and 
antiquarks, makes up a total number of strings equal to 2kn. 



3 



The hadronic multiplicities of the strings in (|]) are obtained from a convolution 
of momentum distribution function and fragmentation functions (eqs. (3.1) to (3.4) of 
0). The former are given by a product of Regge propagators times a delta function 
of energy conservation: 



2kn 

Pkn(Xl, X 2kn ] X 2 , X 3 ,..., X 2 kn-l) = CkriX^^ xf^X^ 1 ■ ■ ■ X 2kn-l$( l ~ X *) ( 7 ) 

1=1 

Here x\ and x 2kn denote the x-values of the valence quark and diquark, respectively, 
x 2 . . . x 2 kn-i those of the sea quarks and antiquarks and c^n is obtained by normalizing p kn 
to unity. The momentum distribution of a single constituent is obtained by integrating (0) 
over the Xj's of all others. In order to regularize the singularities at ~ all Xi (i ^ 2kn) 
in eq. (0) are replaced by |L0| Xi = {x\ + 4/z 2 /s) 1 / 2 with jj? = 0.1 GeV 2 . This introduces 



a free parameter \x in the model. The diquarks and quark fragmentation functions into 
charged particles were determined from e + e~ and/or pp data. They are given by fLOl 



zD qq (z) = 1.12(1 -zf (8) 

, v 1.12 1.3(1 -z) 2 + 0.05 

zD = zD g (z) = 

qv 9sK ' 1.35 1 - 0.5^ 

Finally, we have taken for the threshold of qq — q strings So = (m p + m^) 2 and for the one 
of q — q strings So = (2m^) 2 , with = 0.33 GeV. The value of k can be determined in 
a generalized eikonal model. Alternatively, we can fix it in such a way that the plateau 
in pp agrees with the experimental value of dN ch / dy(y* = 0) for non-diffractive events. 
These two determinations are consistent with each other |7|. |8j . 

Note that at asymptotic energies, when the plateau height of all strings coincides 
(i.e. N qq ~~ qv (y = 0) ~ N q ~ qs (y* = 0)) one recovers the AGK cancellation, namely 



^A(b,y* = 0) ~ n(b)2k N q ^(y* = 0) 
ay 

„ A* T AA (t) a „ D dN^ {v , A> T AA ( b ) = 

°AA{b) pp dy cr AA (b) dy 

Already at RHIC energies, where N Qs ~ qs (y* ~ 0) is substantially smaller than N qq ~ q (y* ~ 
0), the last term of (^) turns out to dominate. However, energy conservation constraints 
produce a decrease of N qs ~ q "(y* = 0) with increasing centrality and, thus, the scaling in 
the number of binary collisions is not reached. Since the pp rapidity plateau dN pp D /dy = 



(1/ (jp p D )d<Tp p D I dy refers to non-diffractive events, for consistency we have to use also the 
non-diffractive cross-section <jp D (rather than the inelastic one) in the determination of 
ua and n, in order to recover eq. @ asymptotically^. 

We can now compute the centrality dependence of hadronic multiplicities. We 
present first the results obtained in the absence of shadowing, at three different energies : 
^fs = 17.3, 130 and 200 GeV. The corresponding non-diffractive cross-sections are JTT], |T2| 
<jnd = 26, 33 and 34 mb, respectively. We take k = 1.4, 2.0 and 2.2, corresponding 
to dN^ D /dy = 1.56, 2.72 and 3.04 [|ll|, |13[. The values of the charged multiplicities 
N^ q (y* = 0) and N^(y* = 0) of the individual strings calculated in DPM, are listed 
in Table 1. The results are shown in Fig. 1 (solid line) and in Figs. 2-3 (dashed lines). 
At SPS we obtain a mild increase of the multiplicity per participant consistent^ with 
the results of the WA98 collaboration || (see Fig. 1). This increase gets stronger with 
increasing energies (Figs. 2 and 3). 

It should be noted that the absolute value of the multiplicity at SPS energies 
determined in DPM has some uncertainty. First, the excess of K + and p over K~ and 
p is not properly taken into account in the above calculation. Second, the value of the 
multiplicity in the q s — q s strings is affected by the value of its threshold. The latter has 
much less effect at -^i =130 and 200 GeV. However, shadowing corrections, which are 
negligibly small at SPS energies, become important at RHIC and have to be taken into 
account. This introduces some uncertainty at RHIC energies (see below). 



As emphasized in ||15|| , shadowing corrections in Gribov theory are universal, i.e. 
they apply both to soft and hard processes. They are closely related to the size of 
diffractive production and, thus, are controlled by triple Pomeron diagrams |15], [16|]. In 



the recent papers [17] it was shown, that, when unitarity corrections are consistently taken 
into account, one can describe hard diffraction at HERA and soft one (photoproduction) 
with the same value of the triple Pomeron coupling. This value agrees with the one in 
ref. [[16|] and will be used in the following calculation. The reduction of the multiplicity 
resulting from shadowing corrections is given by the ratio [T5J 



p fh s Jd 2 s f A (s) f B (b-s) 

Rab{o) = ttt 10 

Tab[o) 



2 Conventionally one uses instead Oinei- In this case one should also use the pp multiplicity for inelastic 
events. Since the latter is about 10 % smaller than the non-diffractive one [jfl], the two effects tend to 
compensate with each other and only a small difference remains in the calculated multiplicities. 

3 The opposite claim was made in jl4| based on a model which is an over-simplified version of DPM. 
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where 



Here the function F is given by the integral of the ratio of the triple Pomeron cross-section 
d 2 a ppp / dydt at t = over the single Pomeron exchange cross-section a p (s) : 



1 d 2 a ppp 

F(s) = 4vr / dy 



(jp(s) dy dt 



C [exp(Ay max ) - exp(Ay min )] (12) 



t=o 



with y = £n(s/M 2 ), where M 2 is the squared mass of the diffractive system. For a particle 
produced at y cm = 0, y max = |M^) ancl Vmin = £n(R A m N /y/3). Using the parameters 
of |T§ we get : C = 0.31 fm 2 . 

The values of the shadowing corrections at each impact parameter R(b) in Au-Au 
collisions at ^J~s = 130 GeV are shown in Table 1. Our results including shadowing are 
shown by the upper lines of the dark bands in Figs. 2 and 3. 

It should be stressed that our calculations refer to dN/dy while the first RHIC 



measurements (18| refer to dN/drj. The latter is, of course, smaller at mid rapidities. 



This difference is negligibly small as SPS where the laboratory pseudo-rapidity variable 
is used. However, at yfs = 130 and 200 GeV, where r\ cm is used instead, their ratio is as 
large as 1.3 ||19|||^. This value is substantially larger than the value of 1.1 quoted in ||. 

Our results for the centrality dependence of dN/drj per participant at s/s = 
130 GeV are shown in Fig. 4 (upper line of the dark band) and compared with the 
PHENIX data. As we see, the centrality dependence is quite well reproduced. However, 
the absolute values are about 15 % higher than the data. 

It should be stressed that the values of R, eq. fllOD are quite large (see Table 2). As 
pointed out in |13|] they have a rather large uncertainty at RHIC energies. An alternative 
calculation of R in ref. |]I3" |, based on a formalism that reproduces the nuclear effects 
in DIS on nuclei, led to values of shadowing about 15 % larger than the ones obtained 
here. Clearly, with these larger values of the shadowing corrections we would obtain a 
quantitative agreement with the PHENIX data. 

Multiplying the values of Ra u Au{b) in Table 2 by a factor 0.85 we obtain the lower 
lines of the dark bands in Figs. 2-4. These bands can be regarded as an estimate of 
4 The increase of < pr > with centrality is predicted to be rather small at RHIC ||(J. Thus, we have 
used the same 1.3 reduction factor for all centralities. 
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the uncertainties on the values of the calculated multiplicities due to uncertainties in the 
shadowing corrections. Note that these uncertainties affect mostly the absolute values of 
the multiplicities - while their centrality dependence is determined quite unambigously 
within our model. 

The multiplicity per participant at SPS energies increases by 1.15 between b = 
10 fm and 6 = 0. The corresponding increase at yfs = 130 GeV and 200 GeV is 1.30 
and 1.31, respectively. As a consequence of this saturation, the rise of the central plateau 
in Au-Au collisions between these two energies is close to the one in pp collisions - an 
interesting prediction of our model. 

We would like to compare our results with the ones obtained in other approaches. 
Many Monte Carlos based on or inspired by DPM or QGSM do contain a term propor- 



tional to the number of binary collisions. However, in other approaches |Ti], [22| such a 
term is associated with minijets. Of course minijets are produced at high energies. They 
have been incorporated in DPM and modify the px dependence of the model | 2q| . How- 
ever, they do not affect the multiplicities, since they play the same role as q s -q s strings 
and the total number of such strings is controlled by unitarity. A comparison with p2| 
shows that, while in this approach the multiplicity is given by a linear combination of n(b) 
and n par t(b), with coefficients independent of b, in our case these coefficients decrease with 
increasing centrality. More important, in our approach these coefficients are calculated 
while in |22] they are fitted to the data at each energy. 

An interesting estimate of the centrality dependence of charged multiplicities in 
high density QCD |24| has been presented in f22fl . Surprisingly the result at y/s = 130 GeV 
is almost identical to the one based on the minijet picture obtained in [22|. However, this 
result relys entirely on the logarithmic dependence of the gluon structure function of the 
nucleon on the saturation scale Q 2 sl used in |f22[ (see last paper of ref. |2^] for a discus- 
sion on this point). It is also interesting to remark that the centrality dependence of the 
multiplicity per participant in the minijet model gets stronger with increasing energy 
- due to an increase of the minijet fraction. On the contrary, in high density QCD the 
effect is the opposite one, namely, the (partonic) multiplicity per participant depends less 
and less on centrality when energy increases^. This is an interesting prediction of the high 
density QCD saturation model. However, in the RHIC energy range, from ^Js = 130 to 
5 This is due to the fact that the geometrical factor p pa rt(b) in Q 2 , which depends strongly on centrality, 
is independent of energy, while the factor x G(x), which depends on impact parameter very mildly, 
increases with s. Therefore £nQ 2 oc £n[x G(x)] at very high energies. 



7 



200 GeV, this effect is negligeably small and the centrality dependence of the multiplicity 
per participant is the same at these two energies. Interestingly, the same result is obtained 
in DPM when shadowing corrections are taken into account. 
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Figure Captions 



Figure 1. The values of dN ch / dy / n part versus n part for PbPb collisions at yfs =17.3 GeV 
in the range —0.5 < y cm < 0.5 computed from eqs. (f|) to (H), compared with the WA98 
data 0| for dN/dr]/n part . The difference between dN/dr] and dN/dy is very small since 
the laboratory pseudo-rapidity is used (see main text). 

Figure 2. The values of dN ch j dy j (0.5n part ) for Au-Au collisions at ^fs = 130 GeV in the 
range —0.35 < y cm < 0.35 computed from eqs. (f|) to (||) (dashed line). The upper line 
in the dark band, is obtained after shadowing corrections computed from eqs. (|Tt})-(|T2|). 
The lower one is obtained with a different determination of the shadowing corrections (see 
main text). 

Figure 3. Same as Fig. 2 for y/s = 200 GeV. 

Figure 4. Same as Fig. 2 for dN ch /dr] c m . The dashed line (results without shadowing) 
has been omitted here. The PHENIX data || are also shown (black circles and shaded 
area) . 



Table Caption : 

Values of N qq ~ q and N q ~ q in eq. (|) at = 130 GeV computed in DPM at 
different values of the impact parameter. The values correspond to charged particles per 
unit rapidity in the range —0.35 < y cm < 0.35. The corresponding shadowing corrections 
Rau Au{b) computed from eqs. ( |TU| - |i"2] ) are given in the third column. 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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b (fm) 






Rau—Au 





0.859 


0.345 


0.656 


2 


0.861 


0.347 


0.657 


4 


0.867 


0.351 


0.664 


6 


0.875 


0.358 


0.681 


8 


0.887 


0.368 


0.712 


10 


0.903 


0.381 


0.763 


12 


0.921 


0.397 


0.843 



Table 1 
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